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ABSTRACT 

The fate of massive cold clmnps, their internal structure and collapse need to be characterised to 
understand the initial conditions for the formation of high-mass stars, stellar systems, and the origin 
of associations and clusters. We explore the onset of star formation in the 75 M© SMMl clump in the 
region ISOSS J18364-0221 using infrared and (sub-) millimetre observations including interferometry. 
This contracting clump has fragmented into two compact cores SMMl North and South of 0.05 pc 
radius, having masses of 15 and 10 Mq, and luminosities of 20 L© and 180 L©. SMMl South harbours 
a source traced at 24 and 70 /zm, drives an energetic molecular outflow, and appears supersonically 
turbulent at the core centre. SMMl North has no infrared counterparts and shows lower levels 
of turbulence, but also drives an outflow. Both outflows appear coUimated and parsec-scale near- 
infrared features probably trace the outflow-powering jets. We derived mass outflow rates of at least 
4 X 10~^MQyr-i and outflow timescales of less than 10'' yr. Our IICN(l-O) modelling for SMMl 
South yielded an infall velocity of 0.14kms^^ and an estimated mass infall rate of 3 x 10~^ M© yr~^. 
Both cores may harbour seeds of intermediate- or high-mass stars. We compare the derived core 
properties with recent simulations of massive core collapse. They are consistent with the very early 
stages dominated by accretion luminosity. 

Subject headings: dust, extinction — ISM: clouds — ISM: individual (ISOSS J18364-0221) — ISM: 
jets and outflows — ISM: kinematics and dynamics — stars: formation 
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1. INTRODUCTION 

1.1. Star formation in massive cold clumps 

The early phases of star formation occur in the dens- 
est and coldest regions within molecular clouds, i.e. 
dense cold cores. Detailed studies of nearby star-forming 
molecular clouds have led to a conception of the evolu- 
tion of individual low-mass stars from prestellar cores to 
prc-main -sequence stars (e.g. Andre fc Montm crle 1994; 
Idi Francesco et al.l 120071 : IWard-Thompson et ahl l2007h . 
The formation of high-mass stars (M > 10 M©) is 
less understood, in particular its earliest stages and 
the link to the origin of associations and clusters 
(jZinnecker fc Yorkd[2007t [Bcuther ct al. 2008). Previ- 
ous studies, mostly aimed towards luminous infrared 
sources, have provided evidence f o r coUimated outflows 
(rece ntly e.g. iZhang et all l2007t iBeuther fc ShepherdI 
120051 and ref. therein) and large-scale rotating structures 
around forming high-mass objects (fZhang 2005, and ref. 
therein). This indicates that they also grow via accre- 
tion through a disk and that the high r adiation pressure 
could be released in outflow cavities (|Krumholz et al.l 
|2005() . In the low-mass regime, accr etion rates are 
found to be of the order 10~^ M© yr~^ (I Whitnev et al.l 
Il997( ). Much higher values of 10~^ M© yr"^ a re expected 
during the build-up of high-ma ss stars tMcKee fc TanI 
l2003t iBaneriee fc Pudritd I2007D and corresponding in- 
creased outflow en ergetics have b een observed (e.g. 
see compilation by IWu et all |2004 and ref. therein). 
Recently, simulations of protostars accreting at such 
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high rates predicted a bloating by several tens of 
solar radii or more and the accompanying high lu- 
minosity at low surface temperature before they 
reach the main sequence ()Yorke fc Bodenheimeij 120081: 
iHosokawa fc Omukail [200l . This gives further indica- 
tion that the formation process differs from a scaled-up 
version of the low-mass case. 

It is presumed that the onset of high-mass star forma- 
tion occurs in cores (diameters of O.lpc and less) within 
massive cold clumps (diameters of around 0.5 pc). Sur- 
veys in the far-infrared and at millimetre wavelengths 
are well suited to search for such objects, but so far 
the earliest, i.e. preste llar, stages have eluded discov- 
ery. iMotte et al.l (|2007( ) derive statistical lifetimes of less 
than 10'^ yr for the former, and of lO'' yr for massive dense 
cores that are forming stars but exhibit low infrared emis- 
sion. Interferometric observations provide the essential 
spatial resolution to study the substructure of clumps, 
and to infer the properties of embedded cores further 
requires a good coverage of their spectral energy distri- 
bution. In this work we present such a detailed study of 
one clump. In particular, we investigate the amount of 
fragments in the clump and their properties, the aris- 
ing outflow activity, and the indications for collapse. 
The clump is located in the region ISOSS J18364-0221 
that has been identif ied using the ISOPHOT Serendip- 
ity Su rvey (ISOSS, iLemke et all (|1996f l: IStickel et all 
(|200l ) at 170 /xm t o esta blish far-infrared c o lour t em- 
peratures. iKrausd ()2003l ) and iKrause et all (|2004l ) se- 
lected sources containing a large fraction of cold dust, 
and su bsequent studies revealed a population of cold 
clumps ()Krause et al.ll2003t iHennemann et al]|2008l ) and 
a collapsing massive core (jBirkmann et al.ll2007f r 
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1.2. The ISOSS J18364-0221 star-forming region and 
the SMMl clump 

The star-forming region ISOSS J18364-0221 
(R.A. 18^36 "24.7", Peel. -02 ° 21'49" fJ2000]) was initially 
studied by iBirkmann et al.l ()2006D and is located at 
a distance of about 2.2 kpc. Derived from the near- 
infrared extinction over a 15' x 15' field (~100pc^), 
the cloud complex associated with the star-forming 
region comprises a mass M « 3200 M0. In the same 
manner, a lower mass limit M > 460 M0 was found 
for the central region (~lpc^), while the far-infrared 
and submillimetre measurements give a luminosity L 
« 800 Lq, an average dust temperature of about 15 K 
and a mass of OOO^ggg Mq . This region contains two 
clumps detected in the submillimetre continuum named 
SMMl and SMM2 that are separated by about 1.5' 
along the east-west direction. The SMMl clump is 
subject of this publication. Its effective radius is about 
0.2 pc, and from the thermal emission a characteristic 
dust temperature of I6.5I3K and a clump mass of 
MsMMi = 75 ± 30 Mq were found. The m ol ecular 
line observations reported in IBirkmann et al.l (j2006f ) 
show red-shifted self-absorption that is interpreted as 
signature of collapse motions. Furthermore, at least one 
outflow is present. The outflow properties resulting from 
the single-dish observations, i.e. the outflow mass of 
about 18 Mq and the mass outflow rate of 10~^ Mq yr~^, 
are comparable to values derived for outflows from pre- 
sumed high-mass star precursors. These results render 
the SMMl clump a promising object to study the 
early phases of collapse and fragmentation occurring in 
massive cold cores. 

The detailed study of this source we present here made 
use of recently collected infrared and (sub-)millimetre 
data described in the next section. They reveal two star- 
forming cores and molecular outflows (Section 3). In 
Section 4, we discuss the core properties in the context 
of early-stage star formation and present a comparison of 
the observed HCN spectra with simple radiative transfer 
models. Finally, the results are summarised in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Near-infrared observations 

Imaging observations in the J, H, and Ks band to- 
wards ISOSS J18364-0221 have bee n obtained and are 
described in IBirkmann et al.l ()2006t ) . Additional near- 
infrared images in the H2 = 1-0 S(l) line (A = 2.122 yum) 
were taken with the Calar Alto 3.5 m telescope in Octo- 
ber 2005 using the prime-focus wide- field camera Omega- 
2000 (Baumci sTer et al.ll2003l ). Omega-2000 features a 
field of view (FOV) of 15.4 x 15.4arcmin^ with a pixel 
scale of 0.4496" pix~^. The exposures were dithered on 
source to allow for sky subtraction. The data reduction 
was done using IRAF. 

2.2. Mid- and far-infrared observations 

IRAC (|Fazio et al .l [200l_ imaging in a ll four pho- 
tometric bands, MIPS (Ri eke et al.ri2004l ) imaging at 
24 /im and 70 /xm and MIPS spectral energy distribution 
(SE D) mode observ ations were undertaken with Spitzer 
(jWerner et al.l l2004| ) . The basic flux calibrated imag- 
ing data of the Spitzer Science Center (SSC) pipeline 
were used for further data reduction and analysis. Cos- 



metic corrections and astrometric refinement were per- 
forme d with the MOPEX software (|Makovoz fc Marleaul 
|2005[) . and the final images were combined using scripts 
in IRAF. Aperture photometry and PSF fitting was 
done with the aperture corrections given in the IRAC 
data handbook and on the SSC website^. The MIPS 
SED mode calibr ation is based on a spectrum of a Boo 
()Low et al.ll2005D and the measured MIPS 70 /zm fluxes. 
The calibra t ion u ncertainties are a bout 2% (IRAC, 
i Reach et al.l (HOGS')), 4% (MIPS 24. lEngelbracht eFall 
(20071)), and 10% (MIPS 70, iGordon et al.l (|2007D ). 
The resulting photometric accuracy is estimated to 5% 
(IRAC), 10% (MIPS 24), and 20% (MIPS 70 and SED). 

2.3. Submillimetre observations 

The submillimetre continuum ob servations with 
SCUB A at the JCMT are outlined in IBirkmann et all 
(j2006[) . In light of the results from the inter- 
ferometric observations that are described in the 
next section, we reanalysed the j iggle maps. The 
ORAC-DR dJe nness fc Economoul [l999l) and SURF 
penness fc Ligh tfoot 1921) software were used for data 
reduction and the photometric calibration based on 
maps of Uranus acquired shortly before and after 
the observations. Furthe r analysis as described in 
iSandell fc Weintraubl (|2001h used the MIRIAD software 
(jSault et al.lll995D . The deviations of the JCMT beam 
from a single Gaussian have been considered by using 
the Uranus maps to construct symmetric beam models 
and deconvolve the maps of the target regions. The 
derived beam sizes are 8.2" at 450 /im and 14.8" at 
850 /im. The maps were restored with Gaussian beams 
of 8" and 14", respectively, and fluxes as well as decon- 
volved source sizes have been derived by fitting Gaus- 
sian components. The noise levels (Icr) in the restored 
maps are lOOmJybeam"^ at 450 //m and 23mJybeam~^ 
at 850 /im. The photometric accuracy obtained is esti- 
mated to 30% at 450 /im and 20% at 850 /im. For a large 
aperture covering the S MMl clump, the photom etric re- 
sults reproduce those of IBirkmann et al] (|2006f ) within 
the uncertainty ranges. 

2.4. Millimetre observations 

We have carried out millimetre observations using 
the IRAM 30m and Plateau de Bure Interferometer 
(PdBI). The molecular lines C0(2-l), HCO+(1-0) and 
IICN(l-O) have been observed together with the con- 
tinuum at 1.3 mm and 3.4 mm. The line frequencies 
are 230.538GHz for C0(2-l) and 89.188526GHz for 
HCO+(1-0). The HCN(l-O) transition includes three 
hyperfine components within 4 MHz at 88.630416 GHz 
(Fi^i), 88.631847GHz (Fa^i), and 88.633936 GHz 
(Fp^i)'^. The PdBI configurations C and D were utilised; 
D was observed in September 2006 and C in April 
2007 with the new generation facility receivers. Spec- 
tral resolutions of 40kHz (3mm lines), 160kHz (C0(2- 
1)), 1.25MHz (3.4mm continuum, 3 x 160MHz band- 
width), and 2.5 MHz (1.3 mm continuum, 2 x 320 MHz 
bandwidth) were used. Phase calibrators were 1749- 
096 and 1741-038, additional amplitude calibrators were 

http : //ssc . spitzer . caltech . edu/mips/ape rcorr/ 1 
NIST Recommended Rest Fre quencies for Observed Interstellar 
Molecular Microwave Transitions: |http: //pliysics.nist.gov/restfreq| 
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MWC349 and 3C273. Corresponding short-spacing ob- 
servations were accomplished in February 2007 at the 
IRAM 30m as On-the-fly maps using the single pixel het- 
erodyne receivers with the VESPA correlator and spec- 
tral resolutions of 80 kHz (1mm lines, 80 MHz band- 
width) and 20kHz (3 mm hues, 40 MHz bandwidth). The 
data were reduced and calibrated with the GILDAS'' 
software. GILD AS was also used to combine the short- 
spacing and interferometric data for the lines. In the case 
of CO, we chose to convert the measured visibilities to 
maps using a weighting scheme that achieves the highest 
spatial resolution (synthesised beam size of 1.85" x 1.33") 
at the expense of sensitivity for extended emission. 
Therefore the combined CO map does not recover the 
complete flux. Compared to the single-dish spectra, the 
line wing fluxes in the combined map are lower by a factor 
of 5 ± 1. For HCN, the comparison of the extracted spec- 
tra shows that the flux measured in the single-dish data is 
reproduced in the combined map. The synthesised con- 
tinuum beam sizes are 1.9" x 1.0" (PA 11.6°) at 1.3 mm 
and 4.6" x 3.2" (PA 6.1°) at 3.4mm. We reach noise lev- 
els {la) in the cleaned maps of 0.6 mJy beam~^ (1.3mm 
continuum), 0.13 mjy bearn"^ (3.4mm continuum), 
0.5Jybeam-ikms-i (C0(2-l)), 0.07 Jy beam"! kms'^ 
(HCO+(1-0)), and 0.1 Jy beam"! kms'^ (HCN(l-O)). 
Assuming a dust temperature of 20 K, we are sensitive 
to total masses of about 0.2 Mq (3cr in the 1.3 mm con- 
tinuum). The noise level (Icr) in the individual HCN(1- 
0) spectra used for the modelling is O.IK. The 1.2mm 
continuum was observed with the MAMBO-2 bolometer 
array in March 2007. The MOPSIC software was used 
for data reduction and the noise level (Icr) in the result- 
ing map is lOmJybeam^^. The photometric accuracy is 
estimated to 20%. 

3. RESULTS 
3.1. Multiwavelength maps of SMMl 

In Figs. [1] and [2] we show multiwavelength maps ob- 
served towards ISOSS J18364-0221 SMMl. In the inter- 
ferometric millimetre continuum observations SMMl is 
resolved into two cores SMMl North and SMMl South 
(Fig. [2 lower panels), separated by 9.5" (ca. 21000 AU). 
In the near-infrared (Fig. [T] upper left) , noticeable ex- 
tinction and reddening is present in the surrounding of 
SMMl. There are no near- infrared sources detected to- 
wards the (sub-)millimetre peaks. The mid-infrared mea- 
sured with IRAC (Fig. [T] upper right) is dominated by 
extended emission in the 8 /im band along the outer rim 
of SMMl. Towards the centre of SMMl, a filamentary 
structure is observed in absorption, and no stellar ob- 
jects are associated. At 24 /im (Fig. [1] lower left) a point 
source at SMMl South shows up. At this wavelength, 
the extended emission around SMMl is also visible, and 
in the north-west a dip remains, but the absorption fea- 
ture is superposed by the PSF of the southern source. 
No obvious 24 /im source towards SMMl North is de- 
tected. Similarly, SMMl South emits at 70 /im (Fig. [T] 
lower right). The wide PSF overlays possible extended 
emission at the rim of SMMl, and no obvious emission 
at SMMl North is detected. 

In the interferometric 1.3 mm contimmm map (Fig. [5] 
top row), both SMMl North and South appear slightly 

^ |http : //www ■ iram . f r/IRAMFR/GILDAS | 



extended to the north, but the beam sidelobes may af- 
fect the morphology. At 3.4 mm they appear unresolved 
(Fig. bottom row), and none of the extensions are 
traced. The contours in Fig. [2] derived from the emis- 
sion in the line wings of the C0(2-l) and HCO+(1-0) 
transitions reveal two molecular outflows. SMMl South 
constitutes the origin of a north-east-to-south-west (PA 
ca. 50°) outflow that also gives rise to the mid- infrared 
features. Towards SMMl North an outflow in east-to- 
west direction (PA ca. -80°) is found. Its red lobe blends 
into the red lobe of the SMMl South outflow. Most 
prominent in C0(2-l) is the SMMl South blue outflow 
lobe that appears coUimated. In HCO+(1-0), the blue 
lobe of the SMMl North outflow is strongest and it also 
shows a rather low outflow opening angle, but in both 
cases the outflows are traced only close to the cores and 
we could not derive an accurate quantitative measure. 

In Fig. [3] we show the near-infrared H2 1-0 S(l) 
(A = 2.122 /im) line emission map. This map has been 
derived by subtracting the scaled K band image from 
the narrow-band image so that stars cancelled out. Stars 
that were not properly subtracted have been masked af- 
terwards for clarity. As shown in the inlay of Fig. [31 
towards the blue and red lobes of the molecular out- 
flow driven by SMMl South there is patchy near-infrared 
H2 emission. Such features trace protostellar jets, and 
the emission proba bly arises f rom collisional ly excited 
H2 in shocked gas (lEliasI 119801: [Reipurth fc Banvll200lL 
and ref. therein). This is supported by faint filamen- 
tary emission in the 4.5 and 5.8 /im bands we detected 
lateral to SMMl South along north-east-to-south-west 
and in the vicinity of SMMl North (yellowish features 
in Fi g. [1] top right). These bands contain several H 2 
lines (|Noriega-Crespo et al.|[200l iSmith fc Rosenll2005D . 
Only a weak near-infrared feature is detected close to 
the blue lobe of the SMMl North outflow. On the large- 
scale map we found additional filaments of H2 emission. 
Though they are not exactly aligned with the molecu- 
lar outflows, they can be assigned to individual lobes 
which would make it necessary that the proposed jets 
show bending resulting in S-shapes. The outermost fea- 
tures in the North-East and South- West are located at a 
projected distance of about 6.7 pc from each other, and 
the South- West feature is at a projected distance of about 
4.6 pc from SMMl South. Among the sources detected 
in the vicinity of the features, we do not find other can- 
didates for jet-driving young stellar objects. 

3.2. SEDs and core sizes 

In the near-infrared and at 3.6 /im no emission is coin- 
ciding with SMMl North or South and we derived upper 
flux limits from our maps. Close to SMMl North, the 
emission in the 4.5 and 5.8 /im IRAC bands were mea- 
sured on levels of 71 and 155 /iJy. In the 4.5 /xm band 
there is extended emission overlapping with the SMMl 
South position, where the flux is about 1.2 mJy. At 8 /im 
no counterparts corresponding to SMMl North or South 
are detected. 

The source towards SMMl South clearly dominates the 
emission at 24 and 70 /im. Due to the separation of about 
9" between SMMl North and South, the derivation of 
fluxes at these wavelengths is impaired towards SMMl 
North (PSF FWHM of about 6" and 17"). We accom- 
plished PSF subtraction to remove the emission of SMMl 
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Fig. 1.— Multiwavelength 2' X 2' maps towards ISOSS J18364-0221 SMMl. Upper left: Near-infrared colour composite of J, H, and Ks 
band observations with Omega-2000 and solid contours of the SCUBA 450 ^m map (Scr, lOo", . . . ). The dashed contour shows the region 
used for background subtraction in the 450 and 850 fim maps. Upper right: Mid-infrared colour composite of the 4 Spitzer IRAC channels 
and contours of the SCUBA 850 fira map (5cr, lOir, . . . ). Lower left: Spitzer MIPS 24 fim map and contours of the PdBI 3.4 mm continuum 
observations (Scr, 6a, . . .). Lower right: Spitzer MIPS 70 fim map and contours of the PdBI 3.4mm continuum observations (3(t, 6a, . . .). 
The 3.4 mm contours in the lower panels reveal the two cores SMMl North and South. 



South. At both 24 and 70 /im the residual maps do not 
reveal a second compact source at SMMl North. The 
3(7 upper limits for SMMl North were derived from the 
residual scatter at its position. We ended up with 24 fim 
fluxes of 199 mJy and 15.7 Jy for SMMl South and upper 
flux limits of 5.3 mJy and 177 mJy for SMMl North at 
24 and 70 fj,m respectively. 

To derive submillimetre fluxes for the two internal 
SMMl components we fitted a Gaussian to each peak 
on the 450 map and an extended Gaussian to ac- 
count for the surrounding clump emission. The same 
sizes convolved with the 14" Gaussian were then used as 



fixed parameters to extract 850 /im fluxes. The clump 
background is well fit with a deconvolved FWHM size 
of 68" and total fiuxes of 12.3 Jy (450 //m) and 2.02 Jy 
(850 /im). For SMMl South we got a deconvolved 
FWHM size of 5.9" and fluxes of 4.42 Jy (450 /Ltm) and 
622 mJy (850 /im), and for SMMl North a size of 7.0" 
and fluxes of 2.66 Jy (450 /im) and 476 mJy (850 /xm). 
T he total fluxes w e deriv ed are higher than those given 
in lBirkmann et aT] ()2006f ) because we allowed for a larger 
extent of the background component. As in the 850 /im 
map (Fig. [1]), SMMl North and South blend into a sin- 
gle elongated maximum in the single-dish 1.2 mm map. 
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Fig. 2.— ISOSS J18364-0221 SMMl CO(2-l) and HCO+(1-0) observations. Top row: Observations at 1.3 mm; continuum emission is 
shown in grey-scale, and CO(2-l) line emission in dashed red (velocity range 20. . . 30kms~^; contours at 2, 3,..., 7 Jy beam~^ kms~^) 
and sohd blue (40. . . 50kms~^; contours at 3, 4, . . . , 9 Jybeam"^ kms~^). Besides, the beam pattern in steps of 20% and the integrated 
CO(2-l) spectrum are plotted. Bottom row: Observations at 3.4mm; continuum emission is shown in grey-scale, and HCO+(1-0) line 
emission in dashed red (velocity range 20. . . 30kms-l; contours at 60, 80, ... , 200mjybeam-i kms"!) and solid blue (40. . .SOkms"!; 
contours at 100, 120, . . . , 400mjy beam~^ kms~^). Besides, the beam pattern in steps of 10% and the integrated HCO+(1-0) spectrum 
are plotted. 



Therefore we extracted fluxes in the same way for both 
sources, using the positions from the submillimetre and 
allowing for a pointing offset. The derived fluxes are 
151 mJy for SMMl South and 156 mJy for SMMl North. 

The continuum fluxes for both cores are listed in Ta- 
ble [1] Fig. [3] shows the spectral energy distributions of 
the continuum emission towards the two sources includ- 
ing the jet features in the IRAC bands. The curves fitted 
to the data are described in the next subsection. 

From the interferometric continuum maps we also de- 
rived fluxes for both sources by fitting Gaussian compo- 
nents. In the case of SMMl South, 80mJy (1.3 mm) and 
4.1 mJy (3.4 mm) and FWHM dimensions of 4.4" x 3.6" 
(PA -50°) were measured, corresponding to a core size 
of -9000 AU. For SMMl North the fluxes are 49mJy 
(1.3 mm) and 3.9 mJy (3.4 mm) and FWHM sizes of 
6.3" X 3.1" (PA 34°) were measured, corresponding to 
a more elliptical core of 14000 x 7000 AU. 

We also estimated the bolometric luminosities from the 
two SEDs. SMMl South has a luminosity L « 180 L©, 
and for SMMl North we got L sa 20 L© using the upper 
limit fluxes in the infrared. 

3.3. Dust temperatures and masses 

lOssenkopf fc Hennind ()1994f ) have derived dust opac- 
ities for dense protostellar cores and for our analysis 



we used the model with thin ice mantles and coagula- 
tion at a gas density of uh = 10^ cm""^ (0H5, Ki.3mm= 
0.9cm^g~^). To estimate the temperature uncertainties 
we also applied the model without ice mantles (0H2) and 
the initial, non-coagulated opacities (OHl). In Fig.|3]we 
show the 0H5 fit to the emission of SMMl South using 
two modified Planck components assuming optically thin 
emission. The long-wavelength SED is matched by ther- 
mal dust emission at 22l^K (solid line). From the fit 
we derived a dust mass of MJj°^'^ = 0.12Mq. The fiuxes 
at shorter wavelengths are reproduced by a dust compo- 
nent at 48 K with a dust mass of M*™ = 5 x 10"^ M© 
(dashed line). For the SMMl North source we used 
the upper fiux limit at 70 /zm to derive an upper limit 
dust temperature of 15"*"^ K from a single component fit 
also shown in Fig. HI The corresponding dust mass is 
Md = 0.18 Mq. Assuming a canonical gas-to-dust mass 
ratio of 100, the masses of the cores are ca. 12 M© 
(SMMl South) and 18 M© (SMMl North), with an un- 
certainty of about a factor two (see l4.2p . 

Preserving the dust temperatures and masses, the ex- 
pected continuum fluxes at 3.4 mm are 4.9 mJy (SMMl 
South) and 4.8 mJy (SMMl North), thus the measured 
fluxes trace about 10 M© for SMMl South and about 
15 M© for SMMl North as compact components. The 
differences stem from the filtering in the interferometric 
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Fig. 3. — H2 (A =2.122 fim) emission map towards ISOSS J18364-0221. The dotted lines indicate the proposed connection of the emission 
features to the different outflow lobes. The inlay shows a magnification of the region around SMMl. The positions of the SMMl North 
and South cores are marked with filled triangles, and the CO(2-l) line wing emission of the red (dotted contour, 2.5 Jy beam"'^ kms~^) 
and blue (solid contour, 3.5 Jybeam"^ kms~^) molecular outflow lobes are over-plotted (cf. Fig.|2]top row). 



observations. Using these masses and the core FWHM 
sizes of 9000 AU and 10000 AU, we calculated volume- 
averaged densities (uhs) = 3M/(47r FWHM^) which are 
4 X 10^ cm-3 for SMMl South and 5 x 10^ cm^^ for SMMl 
North. The remaining mass of the SMMl clump is lo- 
cated in a more extended envelope. From the dust tem- 
peratures and peak fluxes in the interferometric 1.3 mm 
continuum map we derived peak column densities of 
2.7 X 10^3 cm"2 (SMMl North) and 2.4 x lO^Scm"^ 
(SMMl South). Table [1] summarises the properties of 
the two detected cores. 

3.4. Properties of the molecular outflows 

From our C0(2-l) map we inferred the properties 
of the outflows. For this purpose it is necessary to 
convert the measured CO emission into molecular hy- 
drogen column densities Njja , and we used the rela- 
tion Nh, = 3 X 10^° cm-^K-^km-^sx /Tmbdv from 
lOsterloh et all (fl99l . The mass was then calculated 
as M= ^mHad^ ^NnaAn, where ^ is the ratio of gas 
mass to hydrogen mass (taken to be 1.36), mna is the 
mass of the hydrogen molecule, and Afl is the solid an- 
gle covered by Nhj • The integrated emission in the CO 
map towards SMMl traces 4.8 Mq, compared to a mass 



of 75 ib30M(7) derived from the dust continuum emis- 
sion (jBirkmann et al.l 120061 ). This discrepancy probably 
stems from the uncertainty of the above relation in com- 
bination with the missing flux in the CO map (see 12. 4p . 
Apparent also in the integrated C0(2-l) spectrum in 
Fig.m the optical depth in the line is signiflcant. Other 
outflow studies have found C0(2-l) optical depths on the 
order of 10 refe rring to a co mparison with ^^CO emission 
as described in IChoi et all ([1993). In the case of SMMl 
South, the fact that near-infrared H2 emission is observed 
towards both molecular outflow lobes indicates that the 
difference in extinction is not very large. This suggests 
a high inclination i of the outflow axis with respect to 
the line-of-sight. Because we cannot further constrain 
the inclination, we assume i = 57.3° for both outflows 
in the following, corresponding to the mean of a ran- 
dom distribution of outflow orientations. In Table [H we 
list the parameters derived as described in H cnning et al.l 
()2000a^ for the different outflow lobes, corresponding to 
line emission above 3cr. Mean outflow velocities were de- 
rived from the mechanical momenta P and the masses 
using <v>=P/M, and halves of the lobe extents in out- 
flow direction were taken as travelled distance to calcu- 
late the dynamical timescales. We interpret the given 
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TABLE 1 

Continuum flux measurements and derived core properties. 



Parameter 


SMMl North 


SMMl South 


24 fira flux (mjy) 


<5.3 


199 


70 fim flux (mJy) 


<177 


15700 


450 /im flux (mJy) 


2660 


4420 


850 /im flux (mJy) 


476 


622 


30m 1.2 mm flux (mJy) 


156 


151 


PdBI 1.3 mm flux (mJy) 


49 


80 


PdBI 3.4 mm flux (mJy) 


3.9 


4.1 



FWHM size (AU) 


14000x7000 (PA 34") 


10000x8000 (PA -50") 


Luminosity (Lq) 


20 


180 


Dust temperature'' (K) 


15 


22 


NP"^"'^b (cm-2) 


2.7 X 10^3 


2.4 X 10^3 


Mass'^ (Mq) 


15 


10 


(riHa) {cm-3) 


5 X 10^ 


4 X 10^ 



^ Derived from the single-dish continuum flux measurements/upper hmits 
between 70 and 1200 /im.*^ Derived from the interferometric 1.3 mm contin- 
uum map and the given dust temperatures.'^ Derived from the interferometric 
3.4 mm continuum map and the given dust temperatures. 
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Fig. 4.— Spectral energy distributions of ISOSS J18364-0221 
SMMl North (black) and South (red, scaled up by 10). Upper 
limits are indicated by triangles. The dots represent the data 
points obtained with SCUBA (850 and 450 ^im) and MIPS (24 fim). 
Histogram-like bars show the error range of the MIPS SED spec- 
trophotometry (53 to 93 fim) for SMMl South. Stars represent 
emission in the IRAC bands (3.6, 4.5, 5.8 and 8.0 fim) associated 
with the core positions that stems from jet features (see text). The 
J, H and Ks bands were observed with Omega-2000. The lines rep- 
resent flts to the SEDs via modiflcd Planck components (see text). 



masses as lower limits, and as described in l2.4[ the true 
masses are probably higher by roughly a factor of five. 
The same applies to the mechanical momenta and the 
kinetic energies. While the mass derivation may be pre- 
cise within factors of 2 to 4 (cf. [Cabrit & Bcrtout 199(3,), 
the dynamical parameters are less certain and can be 
considered as order of magnitude estimates. Besides, the 
blending of the two red outflow lobes hampered our pa- 
rameter derivation, so we regard the values derived for 
the blue lobes as more precise. 



3.5. HCN 



emission 



towards SMMl 



The three HCN(l-O) hyperfine components Fi^i, 
F2^i, and Fq^i lie within 4 MHz and the expected line 
ratios are 3:5:1 in the optically thin case for local ther- 
mal equilibrium (LTE). The HCN(l-O) map integrated 



over all three components (Fig. [S] left) shows a pro- 
nounced maximum towards SMMl South, whereas only 
weak emission is detected at SMMl North. It is possible 
that the minor peaks that are offset from SMMl North 
stem from the molecular outflows, because it has been 
shown recently tha t HCN can be present there as well 
(jZhang et al.ll2007t ). The HCN(l-O) spectra towards the 
centre positions of SMMl South and North are shown in 
the middle panels of Fig. [5] 

In the case of SMMl North, the signal-to-noise ratio 
of HCN(l-O) is rather low (~ 6) and the line shapes 
are uncertain. As indicated in the plot, the positions 
of the peaks are consistent with a systemic velocity 
vlsr — 36.1 kms^^. The observed line ratios (neglect- 
ing the noise) depart from the expected optically thin 
LTE ratios: Fo^i/Fa^i = 0.6 and Fi^i/Fa^i = 0.75 
compared to 0.2 and 0.6 in the LTE case. Such hy- 
perflne "anomalies" have also been observed toward s 
star- forming clouds in the past ("Walmsl ev et al.l ll98^. 
'Gonzalez- Alfonso & Cernicharo (1993) have investigated 
this numerically and found that the velocity structure 
in the emitting cloud core could strongly affect the line 
ratios. In particular, they found that an increased 
Fi_,i/F2— ,1 ratio indicates inward motion in the core. 
Line ratios similar to the SMMl North spectrum were 
predicted for a subset of their models. 

In the case of SMMl South the HCN(l-O) spectrum 
is detected with a higher signal-to-noise ratio of up to 
40. By fitting the line wings with a symmetric profile, 
we get vlsr = 34.8 km s~^. The hyperfine components 
are characterised by complex line shapes and cannot be 
described within the optically thin LTE approximation. 
The main features of the line profiles are strong dips 
near the line centres. In essence, such dips can be in- 
terpreted: 1) as an artifact of the data combination, i.e. 
the short-spacing flux is not recovered when combining 
the interferometric and single-dish data; 2) as a physical 
effect, namely, as a result of the self-absorption of the in- 
ternal radiation in the envelope. To illustrate that this is 
not an artifact of data reduction we show the single-dish 
and interferometric HCN (1-0) spectra (integrated over 
20" X 20") before their combination in the right panels 
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TABLE 2 

Outflow parameters derived from the CO (2-1) map. 



Outflow Parameter 


South Red Lobe 


South Blue Lobe 


North Red Lobe 


North Blue Lobe 


Pro j . Velocities (km s ^ ^ ) 


40. . . 50 


20. . . 30 


40. . . 50 


20. . . 30 


Solid Angle (as^) 


86 


195 


27 


119 


Mass'' (M0) 


0.066 


0.15 


0.023 


0.12 


Mech. Momentum^" (M0kms~^) 


0.85 


-2.3 


0.30 


-1.8 


Kin. Energy*^ (J) 


1.7x10^** 


2.5x10^** 


1.8x10^* 


2.4x10^* 


Dyn. Timescale'' (yr) 


7.7x10^ 


7.0x10^ 


3.8x10^ 


5.8x10^ 


Mass Outflow Rate'' (MQyr'^) 


8.6x10-'' 


2.1x10-5 


6.0x10-** 


2.1x10-5 


Mech. Luminosity'' (Lq) 


1.8 


2.9 


3.9 


3.4 



Masses are lower limits because of missing flux.'' Assuming an inclinatii 



57.3° of the outflow axes. 



of Fig. [HI If the dips in the interferometric spectrum 
are a result of the missing short-spacing flux, then we 
should see prominent emission near the line centres in 
the single-dish data. However, both interferometric and 
single-dish spectra have strong dips that indicate that 
these line features are real. 

4. DISCUSSION 

4.1. The fragmentation and column density of the 
SMMl clump 

Our observations reveal two major fragments towards 
the SMM l clump. From its m ean FWHM extent of 
30000 AU (jBirkmann et al.ll2006f ) we calculate a volume- 
averaged density of uh^ = 9 x 10'* cm-^ for SMMl. The 
dust temperature of 16.5 K gives a Jeans lengt h of about 
17000 AU and a Jeans mass of 0.7 Mq (Stabler fc Fallal 
|2005| ). This indicates that the two cores with masses 
more than one magnitude higher are not the direct re- 
sult of thermal fragmentation, because one would then 
expect a nu mber of cores with about one Jeans mas s 
to form (e.g. iJappsen et al.l [2001 iBonnell et al.ll2006f ). 
The projected distance of the two cores is 21000 AU. It 
is rather large compared to the radius of influence of 
radiative fe edback from pro t ostella r objects. The sim- 
ulations of 'Krumholz e t al.l (|2007l ) show that the lat- 
ter is around 1000 AU in the first 20000 yr of core col- 
lapse. This suggests that the two cores evolve individu- 
ally in terms of radiative fe edback, however a kinematic 
influence is not excluded. iKrumholz fc McKeel (|2008f ) 
proposed a lower column density threshold of Igcm"^ 
(Nh2 = 2 X 10^^ cm-^) for the formation of high-mass 
stars. For the SMMl clump, we derived a peak column 
density of 7 x 10^^ cm"^ from the single-dish 1.2 mm map 
peak flux. It does not reach the proposed limit. However, 
the core peak column densities lie beyond the threshold 
(see Table [T]) . This shows that observations with high 
spatial resolution are required to evaluate the core prop- 
erties. 

4.2. Properties of the two millimetre cores 

The two detected cores SMMl North and South ap- 
pear quite similar at long wavelengths because they are 
of nearly the same size and exhibit comparable contin- 
uum emission. However, the derived core masses depend 
crucially on the assigned dust temperatures, and these 
are determined by the flux levels in the far-infrared. At 
wavelengths below ~100/im the SMMl clump of about 
75 M0 is expected to become optically thick for contin- 
uum emission, and our Planck component fit does not 



take this into account. We therefore neglect the fitted 
warm component for SMMl South in the rest of the dis- 
cussion. Also for the cold components, we may have un- 
derestimated the characteristic dust temperature. How- 
ever, because of its mass the far-infrared continuum opti- 
cal depth of SMMl is expected to be of the order of unity 
only. This means that a significantly higher dust temper- 
ature, corresponding to a lower mass and in consequence 
to a lower optical depth, would not be consistent with 
the far- infrared fluxes. Thus, the resulting uncertainty 
of the masses are about a factor of two, but the general 
averaging along the line-of-sight could introduce larger 
errors. Furthermore, our mass derivation relies on the 
used dust opacities. While our choice of the 0H5 model 
is rather conservative and the results can be compared 
to other studies, the true opacities may well differ. In 
the following, we take the derived masses at face value. 

Both SMMl North and South are compact and ap- 
pear to drive individual outflows. Therefore we sus- 
pect them to form individual stars or stellar systems 
of a few. With about 10 and 15 Mq, SMMl South 
and North represent cores of intermediate mass. From 
the comparison of core and stellar mass distributions, 
there are tentative indications that the star-forming effi- 
ciency on the scales of individu al cores is between 25 and 
50% (see e.g. lKrumhoI9 l2008l and ref. therein). Under 
this assumption, one would expect both cores to form 
intermediate-mass protostars or protostellar systems in 
case none of the surrounding matter enters in the core 
collapse. However, the latter is probable because of the 
large-scale collapse motions that have been observed to- 
wards SMMl. In comparison to samples of low-mass 
dense cores, S MMl North and S outh are denser by a 
factor of 10 (cf. iMotte et al.]l2007t l. We regard the peak 
column density values we derived from the interferomet- 
ric 1.3 mm peak fluxes as lower limits because of the fil- 
tering in the interferometric map. For low-mass dense 
cores, pea k column densities of lO^"^ cni -^ or lower were 
reported (jWard-Thompson et al.l I1999D . so both cores 
lie above. They rather resemble the objects found in 
IRDCs bv iRathborne et al.l (|2007l ). While their sample 
of four clumps (we use a different nomenclature here) 
have masses of 100 M© and greater, the small scale cores 
are several thousand AU in size and their masses are 2 to 
21 Mq with the exception of one hot molecular core can- 
didate. As in the case of SMMl South, 24 jjum sources are 
associated w ith those co res. Similar findings have been 
reported by lBeuther~fc^ S tcinackcr (2007) for the more 
massive core in IRDC 18223-3, in particular the remark- 
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Fig. 5. — ISOSS J18364-0221 SMMl HCN(l-O) observations. Left panel: The 3.4 mm continuum emission is shown in grey-scale, and 
the contours give the integrated HCN(l-O) line emission at 0.1, 0.2,..., l.ljybeam~i kms-i. The middle panels show the HCN(l-O) 
spectra towards SMMl North and South from the combined IRAM 30m and PdBI data. The hyperfine components are marked for 
VLSR = 36.1kms-i (SMMl North) and vlsr = 34.8 kms'^ (SMMl South). In the right panels the IRAM 30m single-dish spectrum and 
the IRAM PdBI spectrum before combination are plotted. There, the hyperfine components are marked for vlsr = 35kms~^. 



able low luminosity. 

The SED of SMMl South at wavelengths below about 
80 fj,m (Fig. 2]) does not resemble the single thermal 
emission component reproducing the submillimetre and 
millimetre fluxes. A similar result was also derived 
in the study of a sample of 12 cl umps detected in 
five o ther ISOSS star-forming regions (jHennemann et all 
^00§) and for one core in the ISOSS J23053-I-5953 region 
((Birkmann ct al. 2007). Warm and hot dust components 
are required to explain the emission down to 24 /im, and 
also emission from very small grains may contribute in 
this wavelength regime. The fact that SMMl South ap- 
pears as point source at 24 and 70 fim shows that this 
core contains a compact region of heated dust. Assum- 
ing half of the PSF FWHM (6" at 24 ^m) as an upper 
limit, the emitting region has a size of less than 7000 AU. 
This is consistent with an embedded young stellar pre- 
cursor that constitutes the driving source of the SMMl 
South outflow. 

SMMl North lacks emission in the infrared. The only 
sign of star formation is the SMMl North outflow, indi- 
cating that this core is also further evolved than a sup- 
posed prestellar stage. The characteristic dust temper- 
ature is constrained to an upper limit, and therefore in 
this case the derived core mass represents only a lower 
limit. 

Regarding the motion of the two cores with respect 
to each other, their systemic velocities derived from the 
HCN spectra differ by about 1.3 km s'^ and SMMl North 
appears to be receding. The estimated uncertainty is 
0.5kms~^. Nevertheless, the value is relatively high 
compared to the velocity dispersions of up to 0.5kms~^ 
measured for the core-to-core motions e.g. in NG C 1333 
and Perseus (| Walsh et al.ll2007HKirk et al.ll2007f ). It cor- 
responds to roughly 0.3AUyr~^, and on a timescale of 
about 10^ yr the cores would cover their projected dis- 
tance. 



4.3. The star-forming process within the cores 

The SEDs of both cores are dominated by the emission 
at long wavelengths arising from cold dust. This emission 
does not disclose much information about the internal 
structure of the cores. The lack of mid-infrared emis- 
sion and the low spatial resolution in the far-infrared, 
when compared to the expected scales of emitting re- 
gions, do also prevent clarification of the properties of 
the embedded source in the case of SMMl South. SMMl 
North is very likely less evolved, but we cannot exclude 
projection effects. More observations are needed to fur- 
ther constrain the star formation process in the detected 
cores. The core emission morphology in the far-infrared 
will be explored with the ESA Herschel mission. In the 
case of SMMl South, high spatial resolution imaging be- 
yond 20 /im combined with high sensitivity, possible with 
the MIRI instrument on JWST, will severely constrain 
the properties of the embedded protostellar object. In 
the following, we discuss implications given by the de- 
rived outflow properties and our modelling of the HCN 
spectra. 

4.3.1. Outflow activity 

Collimated outflows have been observed for low-mass 
cores and, more recently, also for young high-mass 
objects, supporting the disk accretion scenario (see 
iBeuther fc ShepherdI l2005l . and ref. therein). In the 
former high degree of collimation is linked to 

early stages and a widening of th e outflow opening angl e 
follows during further evolution (|Arce fc Sargentll2006( l. 
According to some observational evidence, this holds for 
high-mass sources as well and may provide a basis for 
their evolutionary classification. 

The fact that we observed rather collimated outflows 
close to their proposed origins SMMl North and South 
is consistent with the idea that the embedded driving 
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sources have formed and launched the outflows recently. 
This is supported by the estimated dynamical timescales 
of less than 10'* yr. 

The outflow energetics change by orders of magni- 
tude from low- mass cases to luminous high-m ass driv- 
ing sources (e.g. iWu et al.l[2f)nl iHenning et al., .2000bl ) . 
and the span of outflow kinetic energies that have been 
derived is 10^ * J < E ^n < lO^'^J. For the sample of 
iBeuther et"all (|2002bD with L > lO^Lo, they got values 
on the order of Eki„ « 10^^ J. With Eki„ « 5 x 10^^ J, 
the outflows of the SMMl cores reach the same order of 
magnitude, in particular because the adopted masses are 
lower limits. They exceed those of low-mass cores. 

So far we have only considered the molecular outflows 
close to the cores. The traced molecular gas has prob- 
ably been entrained by coUimated jets. If we assume 
that the H2 emission features at larger distances stem 
from these jets, which is supported by the rough align- 
ment with the molecular outflows, we can derive a sec- 
ond timescale for the outflow activity of the cores. How- 
ever, we cannot constrain the jet velocities from observa- 
tions. IBeuther et al.l (|2002b[ ) used a ratio of jet velocity 
to molecular outflow velocity of 20, and for the SMMl 
cores the resulting jet velocities are about 300kms~^. 
This is in the range of velocities that have been ob- 
served for Herbig-Haro fl ows from low-luminosity sources 
(jReipurth fc Ballvll2001f) . From the measured offsets of 
the outermost H2 feature in the south-west we got a 
timescale of approximately 1.8 x 10'' yr for the SMMl 
South jet. This is above the molecular outflow timescale 
of approximately 7 x 10^ yr. For the SMMl North jet, 
the outermost feature in the west gives a timescale of 
about 7 X 10^ yr, which agrees with the molecular out- 
flow timescale. Accounting for the overall uncertainties, 
these values are consistent with the proposed relation of 
jets and molecular flows, which of course does not ex- 
clude other scenarios. 

Parsec-scale jets and outflows have been found towards 
many young low-mass stellar objects, mostly traced by 
optical or near-infrared line emission (jReipurth fc Balivl 
I2OOI). The more energetic outflows in the high-mass 
regime arc therefore expected to also extend to these 
sizes, although few observations were reported to date 
(e.g. lBallvll2008l : IBeuther et ani2002a[ ). The presumed 
jets from SMMl support this idea. The fllamentary and 
patchy structure observed in H2 can be caused by dif- 
ferent factors. Besides the varying line-of-sight extinc- 
tion, the local H2 abundance may play an important role 
in determining where the shock-excited emission arises, 
and the kinematics of internal shocks is presumably in- 
fluenced by the penetrated medium. An interesting pos- 
sibility is the connection to the mass ejection history of 
the driving source, which is linked to the accretion his- 
tory for disk-driven outflows. The several features we 
detected indicate a varying mass-loss rate for the SMMl 
cores. Such a burst mode of acc retion may result from 
disk instability driven by infall ( Voro bvov fc Bas"ul l2006) . 

The S-shape that we observe for the presumed SMMl 
South jet can be interpreted as a precession of the outflow 
axis. The outermost feature we found in the south-west 
approximately lies on the projected outflow axis with PA 
« 50° we derived from the molecular outflow lobes close 
to the core, while the features in between are offset to 
lower PA. This means the core rotation axis is not along 



the line-of-sight, but rather has a high inclination. 

4.3.2. Modelling of the HCN emission 

HCN appears to be more resistant to freeze-out in cold 
dense cores com pared to CO-related molecules as HCO+ 
(jRed man et al.l j2008). Assuming that the HCN emission 
traces the dense gas in the cores, we investigate in the fol- 
lowing if the observed HCN(l-O) spectra can be, in prin- 
ciple, reproduced by simple models for the core. First, 
we considered the "one-layer" model where the core is 
spherically symmetric and homogeneous. We assumed 
that the observed spectrum is obtained towards the cen- 
tre of the core as the antenna beam size is little smaller 
than the spatial extent of the core. The parameters of 
the model are: the hydrogen density upja , the tempera- 
ture Tkin, the molecular column density Nhcn, the tur- 
bulent velocity vturb, and the regular velocity Vj-ad which 
characterises the radial expansion or contraction of the 
core. Note that we do not specify the radius of the core 
as well as the actual abundance of HCN since the emer- 
gen t spectrum for such a core depends on their product 
fcf. iPavlvuchenkov et al.ll2008f) . Therefore, in the frame 
of this model we cannot independently constrain the rel- 
ative HCN abundance, the radius of the core, or its mass. 

Having specified the above parameters, we gener- 
ated the model core and performed a line radiative 
transfer (LRT) simulation with the non-LTE code of 
IPavlvuchenkov fc Shustov (I2004D using molecular line 
data from ISchoier et al.l "( 20051 ). A synthetic spectrum 
of HCN(l-O) was calculated from the LRT simulation. 
In principle, the systemic velocity vlsr is an additional 
parameter. We have assumed the values of 36.1 km s~* 
(SMMl North) and 34.8 km s"! (SMMl South) derived 
from the positions of the spectral features. The consis- 
tency between calculated and observed spectra was eval- 
uated with a x^-criterion and took into account all the 
velocity channels of the combined HCN (1-0) spectrum. 

To search for the best set of the mo del parameters, w ^ 
used Powell's minimisation algorithm (jPress 
We had to specify parameter ranges for the search for the 
best-fit model. In our simulations we used the ranges 
given in the Table|3|to derive the best-fit parameter sets. 
The ranges resulted from several runs where we specified 
parts of the adjacent parameter space and did not find 
good reproductions of the observed spectra. Though un- 
likely, we cannot exclude that additional compatible sets 
of parameters exist. In order to assure that the min- 
imisation routine was not trapped in a local minimum, 
we repeated the calculations starting with several sets of 
initial parameters. The best-fit spectra obtained for the 
one-layer model with the corresponding parameters are 
shown in Fig. |6|in the top row. 

In the case of SMMl North, the observed line ra- 
tios were reproduced fairly well. We found a degener- 
acy between the hydrogen density, the kinetic temper- 
ature and the HCN column density and therefore we 
set UHa — lO'^cm^'^. This resulted in Tkin = lOK and 
Nhcn = 7 x 10^^ cm^^. Furthermore, given the assump- 
tion that the observed HCN(l-O) components are indeed 
single-peaked and relatively narrow, the effect of the reg- 
ular and turbulent velocities cannot be distinguished. 
Therefore we set the regular velocity to zero. The derived 
value of the turbulent ve locity is 0.4 km s~^. Note that in 
contrast to the study of iGonzalez- Alfonso fc Cernicharol 
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we did not have to introduce any regular veloc- 
ity to explain the "anomalous" line ratios. The obtained 
ratios naturally appeared as a result of the non-LTE ex- 
citation in the subcritical density. 

In the case of SMMl South, the one-layer model could 
not reproduce the maximal intensities and the strength 
of the self-absorption dips at the same time. We show 
one representative spectrum from the one-layer model in 
Fig. El 

A natural way to reproduce the observed spectrum 
of SMMl South is to include a low-density envelope 
in the model which should lead to self-absorption. We 
checked this by constructing a "two-layer" model with 
5-1-5 parameters, and chose a configuration where a static 
core (v™Jj° = 0) is surrounded by a spherically sym- 
metric, homogeneous envelope. The increased number 
of parameters makes it difficult to identify a best set 
through the minimisation routine. Therefore we fixed 
some of the parameters based on the following argu- 
ments. First, we assumed that the hydrogen density in 
the core is higher while in the envelope it is lower than 
the critical density for HCN(l-O). Given this assump- 
tion, we expect the formation of strong self-absorption 
dips if the HCN column density in the envelope is high 
enough. We set the densities to n^™ = 10^ cm~^ and 
n™^'^ = lO^cm"^. We also fixed the temperature in the 
envelope = 10 K assuming that it is not heated 

by the inner (most probably warmer) part of the core, 
for which we chose Tj^°j^° = 30 K. Changes in the core 
temperature do not affect the derived model spectra sig- 
nificantly. Thus, we varied the five parameters N^^^, 
'^tSI-b> ^?fcN' ■^turb' "^rad' ^'^^ Constrained them using the 
minimisation routine. 

The best-fit spectrum obtained from the two-layer 
model is presented in Fig. [6] (bottom). The two-layer 
model could reproduce the intensity, the strong self- 
absorption dips and asymmetry of the observed spec- 
trum quite well. The line wings may stem from out- 
flowing gas, or also be due to a more complex velocity 
profile not incorporated in the model. The important 
outcome of this model can be summarised as follows: 
First, the strong self-absorption in the line profiles can 
be reproduced by the high-density core together with the 
low-density envelope. Second, the asymmetry of the line 
profiles can be explained by the infalling envelope (we 
got v™j = — 0.14kms~^), given the relatively high tur- 
bulent velocities in the core {v^°l^ — 1.7kms^^) and in 
the envelope (vj?^^,^ = 0.5kms~^). However, this config- 
uration is not the only way to fit the line asymmetries. 
In particular, we were able to reproduce the line asym- 
metries by the combination of an expanding core and a 
static envelope because of the similar relative velocities. 
We consider this scenario implausible, though, because 
of the short outflow timescales. 

Finally, we stress that the considered simple models 
did not reliably constrain the other parameters of the 
cores, but they illustrated that the observed spectra can 
indeed be explained theoretically. In order to get more 
reliable information about the core properties one has 
to consider a set of molecular lines and transitions to- 
gether with more physical models of the source. Unfor- 
tunately, the available data, in particular the HCO''"(1-0) 
measurements, are not sufficient to enable this. Indeed, 



if we used HCO+(1-0) in addition to HCN(l-O) in our 
fltting routine, we would have to incorporate two more 
parameters for the 1100+ abundance (in the core and 
in the envelope), while IICO''"(1-0) alone will not allow 
us to further constrain the excitation conditions. For 
the latter, at least a few more transitions of HCO+ or 
its isotopomers were necessary. On the other hand, a 
more physical model of the core (e.g. with smooth dis- 
tributions) has more parameters to fit which makes such 
a model even more uncertain. Thus, the usage of more 
sophisticated models is required when detailed spectral 
maps of the sources are available. 

4.3.3. Collapse of SMMl South 

The modelling of the HCN emission showed that the 
spectrum of SMMl South can be explained with a col- 
lapse of the core. From the resulting velocity we cal- 
culated estimates for the mass infall rate using dM/dt 
« M/t = Mvi„/r. For SMMl South, Vi„ = CUkms"! 
is the regular velocity of the envelope in the two-layer 
model. Together with the mass and radius from Table [T] 
this gave an infall rate of 3 x 10~^ M© yr~^. This is con- 
sistent with the mass outflow rates, but is only a crude 
estimate. 

4.3.4. Internal turbulence 

Additional parameters we obtained from the HCN 
emission are the internal turbulent velocities of the cores 
that reproduced the observed line widths. However, in 
both cases rotational motions, more complex velocity 
profiles, and also outflows may contribute to the line 
widths. Such effects were not considered in our mod- 
els. In the case of SMMl South, the fitted model shows 
a very turbulent core and less turbulence in the enve- 
lope. When compared to the isothermal sound speed of 
around 0.2kms~^ for cold cores (T « 20 K), the core tur- 
bulence of 1.7kms^^ is highly supersonic. The turbulent 
velocity in the envelope of 0.5kms~^ is on the scale of 
the sound speed but still in the supersonic regime. The 
SMMl North HCN spectrum indicates a more quiescent 
state of this core. The linewidths are reproduced with a 
turbulent velocity of 0.4kms~^. 

5. SUMMARY AND CONCLUSIONS 

We have investigated the onset of star formation in 
the massive cold clump ISOSS J18364-0221 SMMl. The 
derived results regarding the clump fragmentation, the 
core and outflow properties, and the collapse indications 
are summarises as follows: 

1. Two compact, embedded cores are the only objects 
we identified. Thus, this clump shows little frag- 
mentation, but seems to rather produce a low num- 
ber of compact objects. 

2. The cores dubbed SMMl North and South are of 
about 10000 and 9000 AU radius. The dust tem- 
peratures are 15 and 22 K and they have masses 
of about 15 and 10 Mq, respectively. SMMl South 
harbours an infrared source not detected at 8, but 
at 24 and 70 ^m. The luminosity of SMMl South is 
about 180 L© and a molecular outfiow is detected 
in its vicinity. SMMl North lacks infrared emis- 
sion, but a second molecular outflow is emerging 
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TABLE 3 

Parameter ranges used in the HCN(l-O) 

EMISSION modelling FOR THE ONE-LAYER AND 
TWO-LAYER MODELS. 



One-layer Model 
Parameter SMMl North SMMl South 



iiHs (cm 3) 


10^ 


10^... 10^ 


Tkin (K) 


10 ... 50 


10... 50 


Nhcn (cm 2) 


10^2 ... lO^s 


10^3 . . . 10^5 


vtm-b (kms-1) 


0...1.0 


0. . . 1.0 


Vrad (kmS-1) 





-0.25... 


Two-layer 


Model for SMMl South 


Parameter 


Core 


Envelope 


(cm"3-) 


10^^ 


10^ 


Tkin (K) 


30 


10 


Nhcn (cm 2) 


10^3 . . . 10^5 


10^3 . . . 10^5 


vturb (kms~i) 


0. . .3.0 


0. . . 1.0 


Vrad (kms^i) 





-1...1 



Model illustration SMMl North SMMl South 




20 30 40 50 

Velocity (km/s) 



Fig. 6. — Model illustration, observed and modelled HCN(l-O) spectra towards the two SMMl cores. In the top row, the one-layer model 
is illustrated and the resulting best fit spectra (red curves) are shown with the corresponding model parameters. The bottom row shows 
the two-layer model and the resulting best fit of the SMMl South spectrum (red curve) as well as the corresponding model parameters. 
Parameters that have been fixed are marked with a preceding asterisk (*). 
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from this core. The luminosity of SMMl North is 
approximately 20 . Both outflows appear coUi- 
mated and we derived lower outflow mass limits of 
about 0.2 and 0.3 Mq of presumably entrained gas, 
and mass outflow rates of about 4 x lO^'^ Mq yr~^. 
In Fig. [7] we show a schematic illustration of these 
results. 

3. We detected filaments of shock-excited H2 emis- 
sion at 2.122 /im towards the lobes of the SMMl 
South molecular outflow, and also at distances up 
to 4.6 pc from SMMl. Each filament is roughly 
aligned with one of the outflow axes. They are 
possibly tracing the outflow-generating jets, and 
the dynamical timescales we derived are consistent 
with those of the molecular outflows considering 
the uncertain jet velocities. For both molecular 
outflows, ages of less than 10^ yr were found. 

4. SMMl South is bright in HCN(l-O) and our mod- 
elling results support that the core is collaps- 
ing. We got an infall velocity of 0.14 km s~^. 
This results in a mass infall rate estimate of 3 x 
lO~^M0yr~^. For both cores the spectra can be 
interpreted with turbulent internal motions, and 
in the case of SMMl South the innermost part ap- 
pears highly supersonic. 

Both SMMl North and South are more massive than 
typical low-mass dense cores and embedded in the pre- 
sumably contracting 75 Mq SMMl clump. This suggests 
that the two cores harbour protostellar seeds that may 
become intermediate- to high-mass stars. The presence 
of a 24 fim source, a very turbulent central region and 
the jet features at large distances to the core support 
that SMMl South is more evolved than SMMl North. 
We interpret it as having an embedded young protostar 
that constitutes the driving source of the outflow. The 
outflows as well as the infall indicate that it is accret- 
ing. The associated outflow indicates that SMMl North 
is also a star-forming core, and the outflow energetics 
are similar. However, the forming object remains unde- 
tected, and the relatively low level of turbulence implies 
that it has only marginally affected the surrounding core. 
The presumed difference in evolution of the forming ob- 
jects contrast with the similar sizes, densities and outflow 
properties of the cores. These findings suggest a rapid 
evolution of the lumin osity as predicted for high accre- 
tion rates (|Yorkdl2008D . 

In comparison to t he core collapse simulations of 
iKrumholz et al.l ()2007[ ) we find that the luminosity of 
SMMl South would be reached at a stage where the em- 
bedded protostellar mass is around 0.3 Mq. Accretion 
onto the protostar is the dominant luminosity source, and 
they find an accretion rate on the order of lO"'^ Mq yr~^ 
from directly after the formation of the protostar until 
this stage. The resulting timescale of about 3000 yr is 
consistent with the outflow timescales. Also the mass 
outflow rate as well as the infall rate estimate may be 
consistent with the high accretion rate. SMMl North is 
less luminous, and hence the protostellar mass would be 
even lower, but may still be consistent with a high ac- 
cretion rate in their model. However, the formation of 



outflows is not incorporated therein 



ISOSS J18364-0221 SMM1 
24Ltm & 3.4mm contours 

SMM1 North 

T = 15K 
M = 15M„ 

peak N„ = 2.7x1 0^cm-= 




Fig. 7. — Schematic illustration of the identified cores, the corre- 
sponding outflows, and their properties. The 3.4 mm contours are 
plotted above the 24 /im emission (cf. Fig.[TJ. 
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